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Abstract. Assuming that two incoming annihilating particles interact by exchanging a generally 
massive attractive vector (or scalar) boson, we find, by taking the non-relativistic limit of the field 
theory ladder diagrams, that the complete annihilation amplitude A is equal to: the convolution of a 
solution of the Schroedinger equation (including the attractive potential) with the Fourier transform 
of the bare (i.e. ignoring the attraction) annihilation amplitude Aq. The main novelty is that Aq 
can be completely arbitrary. For a Coulomb potential we find analytically the enhancement for 
the 1-partial-wave cross-section, e.g. the P wave enhancement 2TT{a/vY {v relative velocity), for a 
Yukawa potential we describe a simple algorithm and give numerical results showing an important 
P wave enhancement with a resonant pattern. 

1 Summary and results 

Here we discuss general results for the Sommerfeld enhancement, which could be an important ingre- 
dient for the interpretation of astrophysical data to signal dark matter annihilation processes [I] , [2] , [3] , [1] , 



We assume that two incoming particles attract each other by exchanging repeatedly a massive 
vector boson, before undergoing some annihilation reaction. We will be interested in the non- 
relativistic limit of this process: in this limit the exchange of a scalar boson would give the same 
result, although to be definite we will continue to refer to a vector boson. 

Also, in general one has to consider a non-abelian process, corresponding to a matrix-valued 
interaction. In the analysis below we assume that that the matrix can be diagonalized and that 
there is a boson of definite mass mediating the attraction for the relevant eigenstates of the two 
particles. We consider therefore such a definite eigenstate. 

Let p be the CM momentum of the incoming particles and m their mass. We call AqIp) the bare 
(i.e. neglecting the vector field attraction) annihilation amplitude and A{p) the complete (including 
the effect of the attraction) amplitude. 

Our general result is a generalization of the enhancement formula of ref. IJ. Our treatment 
focuses on the relation between A{'p) and Ao{'p) for a totally generic Ao{p), rather than taking for 
instance a constant Aq (i.e. its Fourier transform ~ ^(r)). This allows us in particular to consider 
the specific case of partial waves higher than S. 

Our method is also similar to the one of ref.[l , even if the path we follow is not quite the same. 
In a sense it is more direct, since for simplicity we focus on the essential points, leaving off the details 
of the dark matter properties that are not relevant for our generic discussion. 

We begin by writing an integral equation for A(p), whose iterative solution corresponds to the 
sum of the ladder diagrams in which the incoming particles repeadly exchange vector bosons before 
annihilating. Then we take the non-relativistic limit finding an in-homogeneous Schroedinger equa- 
tion relating the Fourier transform of A{p) to the Fourier transform of Aa{p). The solution of this 
equation gives: 
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where i^p-, normalized as j dr(j)*^{r)(j)^(r) = (27r)^(5(p — k), solves the Schroedinger equation 

i-T^d'' -^^-^)(j)p{r)=Q (1.2) 

rrir = m/2 is the reduced mass, the vector potential mass and a = g^/(47r) its strength. 

This formula for the complete amplitude A can be numerically evaluated once the dependence 
in p of Aq{p) is known. 

Furthermore, we consider the case in which the annihilation occurs at a definite partial wave I 
and therefore assuming that for p small Aq ^ p\ i.e. the typical dependence on p of the I partial 
wave. 

We define the enhancement enhi as the factor that multiplies the bare cross-section cto,; to give 
the complete one ai (for the rate it would be the same): 

ai = enhi ■ ao,i (1-3) 

We find that enhi depends, beside Z, on the two dimensionless parameters: 

a:^^, b:^J^ (1.4) 
V rrirV 

where rrir = m/2 and v is the relative CM velocity. 

Analytical computations are possible in the case in which ^ is negligible. We find for the generic 
I partial wave : 

In particular, in the limit a » 1, for the S wave we find the standard Sommerfeld enhancement 
27rf , [9],[1],[2],[3], [5], [6], whereas for the F {I = 1) wave we find the enhancement 27r(f )3. 

In the general case of a Yukawa potential, we derive a formula expressing enhi in terms of the 
asymptotic behavior of the partial wave solution of the homogeneous Schroedinger equation, and we 
indicate a very easily implementable algorithm to evaluate it numerically, which does not require 
any computational skill. 

The case of the S Z = wave has been discussed at length in the literature, and it has been found 
a resonant pattern see in particular refs. [6] and |8^. 

Here we present some numerical results for the P I — I wave, which also shows resonances. In 
this case the enhancement can be of several orders of magnitude also away from the resonance, and 
much more on it. 



2 Derivation of the results 
2.1 The equation for the amphtude 

Define A{p,p'; Po) to be the amplitude for the annihilation process of two x particles 

xiPi) + X{P2) a{p[) + aijp'^) 

The final state a, a can be any (quantum- number compatible) two particle state of the standard 
model. 

Define P = = , p = p' = in the CM Pq = ^p^+m^, p ^ q, po 0. 

A(j},p'; Pq) is the complete amplitude, including the Sommerfeld effect, and we call Aq{p,p'; Pg) 
the "bare" amplitude, that is neglecting the Sommerfeld effect. The variable p' does not play any 
role in the following, we continue to indicate it just for completeness. 
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Here we treat the initial particles x having a mass m as Dirac particles of opposite charge, to be 
definite, and imagine that before annihilating they attract each other by the exchange of a vector 
boson called VB (the exchange of a scalar would give the same non-relativistic potential). We take 
the general case of the VB with mass /i << to. In the results one can put /i = 0. In the non- 
relativistic limit, any kind of x particle, Dirac, Majorana or even scalar, would give the same result. 
As said, in the case of non-abelian interaction, we assume that it has been diagonalized. 

Here we follow closely the Chapter 10 of the book "Quantum Field Theory" by C.Itzykson and 
J.B.Zuber [TUj. 

Since there are two fermions in the initial state we distinguish the Dirac matrices acting on the 
two particles by a sufhx, say 71,72 act on particle 1 and 2 respectively. 

The VB exhange between the two incoming particles gives a factor (vertex-propagator-vertex) 

7i72 - 71-72 ^ 7i72 ^ , 7172^0 7i ■ 72 . 

fc^-M^ fc2^_^2 ^(fc2_^2)(fc2+^2) fc2-/i2'' 

where k — p — p'. The first term in the r.h.s. is the instantaneous "Coulomb like" interaction and 
it will treated non-pcrturbatively, whereas the second term, containing the retarded and magnetic 
effects can eventually be included as a perturbation [10] and it will be ignored in the following 
discussion. 

Also, we will include only the ladder diagrams, that is the iteration of the one VB exchange [TO] , 
The reason for it is that the complete amplitude A(p,p'; Pq) will be determined, in the relevant 
non-relativistic approximation, by the non-perturbative solution of a Schroedinger equation with a 
potential to be determined by the field theory diagrams. A usual strategy for solving the Schroedinger 
equation is, whenever possible, to split the potential into a dominant part which can be exactly 
solved, plus higher order terms giving corrections to be computed, if necessary, by perturbation 
theory, like it is done for the fine and hyperfine corrections to the energy levels of the hydrogen 
atom. 

In our case the general diagrammatic expansion of the amplitude can be seen as the iteration 
of two-particle irreducible sub-diagrams and the non-relativistic potential corresponds to the sum 
of those two-particle irreducible diagrams. The dominant term of the potential, for which we will 
provide the exact solution, comes from the lowest order two-particle irreducible diagram, that is the 
one VB exchange. 

In the general expansion it appears also the iteration of other two-particle irreducible sub- 
diagrams, including radiative corrections and sub-diagrams in which the VB lines cross; as said, 
they correspond to additional contributions to the interaction, which are of higher order in the cou- 
pling constant, i.e. higher than the one-VB potential, and their effect , if necessary, can be computed 
by perturbation theory. This is similar to the treatment of the positronium in ref. [10| , in which the 
Bethe-Salpeter equation in the ladder approximation provides the leading non-perturbative solution, 
with perturbative corrections coming from higher order diagrams. 

Further, we also note that the crossed diagrams are less singular for k^ ^ 0{fi^) than the one-VB 
exchange; since the range of the interaction is determined by the strength of the nearby singularity, 
they will correspond to an interaction which is less effective in attracting the incoming particles 
at large distances, beside giving a small correction to the exact solution solution of the one-VB 
potential, because they represent higher orders in the interaction. 

That being said, A satisfies the following integral equation {p = 7'^p^): 

A{p,p';Po) = Ao{p,p';Po)- (2.1) 

. 2 f d^qdqo 7i°72° (P + g + TO)i {P-q + m)2 , 

^ J (2^)4 (^_g-)2+^,2(p + g)2_^2+,,(p_g)2_^2+,, l^^'^' O) 

Its iterative solution corresponds to the sum of the ladder diagrams in which the incoming particles 
repeadly exchange VB's before annihilating. After solving the equation one has to put p = p, that 
is po = = 0, p = p. 
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We evaluate the integration over qo by closing the contour in the complex plane disregarding the 
possible singularities of A which give sub-leading terms in the non-relativistic limit (see Appendix 
A). In the following v := \v\ denotes the modulus of a 3- vector. 

The poles in the lower plane are located at : qo = to ^ Pq — ie, qo = lu + Pq — ie with 
CO = \/ + frP- ■ We write [10 

(P + g + m)i (P-g + m)2 



(P + q)^ — + ie {P — q)^ — 7ti? + ie 



(2.2) 



90 + Po - + «e go + -TO + - «e qo - Pq ~ ui + le qo — Pq + u; — le 
where A^ = and H = a- q + f3m with /3 = 7", a = /37. 
The residue at qo ~ lu — Po is 

A+r,-^r ^2'(-g) A+(-g) A+(g)A+(-g) ^o_o 

here we take the leading term in the non-relativistic limit, that is the term containing at the 
denominator cj — Pq = \/ q^ + rn? — y/p^ + ra^ (small for m — > 00). 
The residue dX qo = uj + Po is 

here there is no small denominator, and we neglect it as a sub-leading term. 

In the non-relativistic limit 7" 1 and A]^(g)Aj(— (7) ~ 1 (leading term for ni 00). Futher, 
in the nonrelativistic limit, A{uj — Pq, g, p'; Po) ~ v4(0, g, p'; Po). Therefore we can consistently put 
directly po = and consider the equation for A[p,p'\ Pq) = A(0,p,p'; Pq). 

In conclusion, in the non-relativistic approximation we get 

Al~ ' P\ A (~ ' .g^ f rf^g A{q,p'-Po) 
A{p,p;Po)^Ao{p,p;Po) + -pr^ — t^tt wr- (2-5) 

[2TTy J (p _ g)2 _,_ ^2 2(w - Po) 

In the denominator oj—Po we take the leading term for m large compared to the three- momentum: 



where £ = 2(Pq — to^)/(2to) = p^/(2mr) is the total non-relativistic energy (m^ — m/2 is the 
reduced mass). 

The final step is a further redefinition: 

A(^,p';Po) ^{^- £)^s{p,p') (2.7) 
getting the non-relativistic equation 

(^-5)v^.(;?,p')-^o(|,p';Po)+ (1^/ -J^I^M^^^P') (2-8) 

Since = j;^ J df " "'J' , defining a := g^/(47r) and 

^£{r} / #e'^'-^"*£(p,p'), C^o(r) :- / d0e^P--^Ao{0,p' , Po) (2.9) 



we get the in-homogenous equation: 



1 fyp~>^^ 

{- £)*£(f) = C/o(f) (2.10) 

Zirir r 
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2.2 Solution of the equation 

Ea. (|2.10p is formally solved by 

r r dk <t>dr)(t>*Ar') 

(we choose the Feynmann contour prescription: — ie, for our purpose +ie would give the 

same) where 4>^{'r) are a complete set of solution of 



(-77^52 - ^i^— - -^)0r(f) = (2.12) 



normalized such that the completeness relation is J (2*^)3 '/'fc(^)0g(^) = <5(r — r*). 
The final step is to reconstruct Aon^sheuij>TV' \Pa) = ^^i^p^f;^{PiP''i Po) on the mass-shell, that 
fo"^ 2^ ^ ^ = 2^ i-e- ^ Pq -rr? ^ . We get 



sheiiip.p'; Po) = liin ~ S)^£{P,p') (2.13) 



dk [/ rfre-P->,-(rO] [/ dr'^;-(r'){/o(r-)] 



(27r)3 ^ V (27r)3 fc2 - p2 _ 

We use the formula (see Appendix B) 

1 ,~2 ,,Id^e-^^ir) 



(2.)3 1^ ^-^' y^-^^ie - '^P- ^) (2.14) 
Then the scattering amplitude A on shell including the Sommerfeld effect turns out to be 

Aon-sheiiip,p';Po) - J dfrpir) J ^e*f Mo(9,y;Po) (2.15) 

This result can be rewritten in a trasparent way as 

Aon-sheii{p,p';Po) = J dq(j)p{q) < q,-q\M\p' > (2-16) 

where < —q\M\p' >= Ao{q,p') is the matrix element of the annihilation reaction, and 4'p{q} 
is the momentum-space wave function of the incoming pair, which takes into account the mutual 
interaction, normalized such that in absence of interaction A — Aq. (A formula of that kind is 
heuristically presented in the Chapter 5 of the Peskin and Schroeder book [TI] . With our Feynmann 
graph convention the initial state appears to the left rather than to the right). 

Let us see how the above general formula works in two particular cases. 

In the case in which Aq is S-wave dominated and then it is a constant, we have: 

Ao{p,p'; Po) = «o ^ 1 ^e'«"%(g,p'; Po) = S{r)ao (2.17) 

Therefore by eq. (|2.15p 

Aon-sheuip,p'; Po) = <t'fiO)ao (2.18) 
In the case AQ{q,p'; Pq) is P-wave dominated then it is linear in p 

Ao{q,p'; Po) = p • f/ai ^ J ^e'^"Mo(g,p'; Po) = -tdSif) ■ p'ai (2.19) 
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therefore by ea. (|2.15p 

Aon-sheuip,p';Po) = id4>*j;iO) ■ jf ai (2.20) 

To make explicit computations, one has to find (f> around r — 0. In order to do that, it is 
convenient first to speciafize the general formula (|2.15p for the case of a definite partial wave I. 



3 The enhancement formula for generic 1-waves 

The wave function 4>p{r) can be decomposed in partial waves [S]: 

Mr) = ^^J^^*'(2; + l)e'''i?p,;W^/(p-^) (3.1) 

(p, f are unit vectors in the direction of p, r). 

Rp,i{r) is the solution of the partial wave Schroedinger equation: 

normalized such that 

r^drRgj{r)Rpj{r)^Sip~q) (3.3) 



and the completeness is 



f°° 1 

dpRpAr)RpAr') = ^S{r ~ r') 



(3.4) 



(This can be checked to be consistent with the normalization J dr4)p{f)4)j:{r) — {2tt)^6{p~ k) and 
with the completeness / d^k(t)^{r)(l)^{r') — (27r)^(5'^(r — r')). We have taken the convention that 
Rp,i{f) is real, which we can always do. 

For the free case (tPp{r) = e'^'' and i?^ ;(r) — lima-^o Rp,i{r)- 

We use the following identities 0: 

An 

dQkPiik ■ f)Pv{k ■ f') = l^j^^ii'Piif ■ f) (3-5) 

y(2; + l)P,(f-f') =2(5(1 -f-f') 5^(f-r')^^5ir-r') — 5(l-f-f') (3.6) 

Take for the I wave : 
Ai{p,p') = Ai{p,p')Pi{p-p') and similarly A^^q^p') = Ao,i{q,p')Pi{q ■ p'). 

Putting the expansion eq. (|3.1|) in the main formula eq. (|2.15p and using eq. p.Sp to do the angular 
integration, we get the partial wave version of eq. (|2.15p 

Ai{p,p') = - r^drRp^r) qdqRli{r)AoM.P') (3-7) 



P Jo 

The standard dependence on small q is Aoj{q,p') — q'-ao,i{p')- Since [9] 



using the completeness eq. (|3.4l) for the free i?° ; we get 



qdqRlir)q' a,,{p') = H'a/| ^ ' ^ ' ^ ^ W a,,{p') (3.9) 
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In conclusion we get 



Mp,p') = (-)' 



7rl-3---(2/ + l)l 
Tl p 



TT 1 • 3 - • • (2? + 1) 1 d I 



r'^drRpj{r)—5\r) aoAp') 
2 II p i^y^r^M^o aoAp') 

4 The enhancement for the case of a Coulomb potential 

The Coulomb partial wave is [9], with our notation a :^ a/v (v relative CM velocity), 



(3.10) 
(3.11) 



- TT (s2 + a2)e-r(l - ta) ^ ^ ' F{ia + ? + 1, 2/ + 2, 2ipr) (4.1) 



It is easy to evaluate {-^Y Rp i{r)\r=o by using the explicit expression ea. (|4.ip since for r ^ 
only the Z-time derivative of the factor (2pr)' contributes and F{a,b,2ipr)\r=o = 1- 
We get El : 



nLi(s' + a')e'^r(l-ia) 



P^a,o,i{p) 



(4.2) 

Since p^aQ^P') = Ao,i{p,p'), and |r(l - i/pc)\'^ = sinllTra) taking the square modulus we get 
the enhancement formula ai = enhi ■ ctq,; where 



enhi = J|(s2 +0^)6" 



sinh{TTa)l 



|2 



(4.3) 



For large a = a/v, enhi = ^^(f)^'^^ 



5 The enhancement for the case of a Yukawa potential 



Here we take ea. p.lip and insert for i?p,;(r) the solution of ea. p.2p . normahzed as eq. 

We know from textbooks (see for instance W) that this normalization corresponds to the asymp- 
totic behavior 



2 sin(pr - ^ + 5;) 



Let us define x = pr and put Rp,i{r) — Npx^^i{x); the equation for is: 



2(^ + 1)^, 2ae-'"' „ 
+ -^'i + ( + 1)*; = 



X 



X 



where again a :~ a/v, b :~ p,/{mrv) and v = p/riir is the relative velocity. 
Suppose we solve this equation with the initial conditions 

$K0) = 1 ^'i{0) = -a/{l + l) 



(5.1) 



(5.2) 



(5.3) 



(the condition for $J(0) is dictated by the equation for a regular solution). Then the asymptotic 
behavior will be 

x'+'M^),^^^Csm{x-^-^ + Si) (5.4) 



^checking the normalization: for o ^ one gets A;(p,p') = p'ao i(p') 
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In order to agree with the normaUzation of eq. (|5.ip we have to put TV = \j '^ZJ- Substituting in 
eq. p.ll|) we get 

Ai{p,p') = ^ -pao.i = ^ -Ao.i{P,P) (5.5) 

In conclusion, by defining the Sommerfeld enhancement enhi for the I partial wave cross-section (or 
equivalently for the rate) as 

ai = enhi ■ ao,i (5.6) 

we get 

enhi = ( ^ '-f (5.7) 

where C is obtained by looking at the asymptotic behavior eq. (|5.4p of the solution of ea. (|5.2p with 
the initial conditions eq. (|5.3p . enhi depends on the two parameters a and b. It is not necessary to 
determine Si. 

Another equivalent strategy is to put Rpj{r) = Npipi{x)/x; the equation for ipi is: 

^r + (l + -e---^)^. = (5.8) 
X a;^ 

If one solves this equation with the initial conditions corresponding to 

Mx)^_^,^x'+' (5.9) 

then the asymptotic behavior will be 

Mx)x^oo^Csinix~i^+ Si) (5.10) 



with the same C of eq. (|5.4p giving the enhancement as in eq. (|5.7l 



5.1 Computations for / = 1 

In principle it is easy to get C: for instance one can use the NDSolve instruction of Mathematica to 
get the numerical solution of eq. (|5.2p with the initial conditions eq. (|5.3p . or equivalently of eq. (|5.8p 
with initial conditions eq. (|5.9p . In order to find C one takes Fi{x) = x'-^^^i{x) or else F{x) = ip{x), 
and one plots Fi{x)^ + Fi{x — 7r/2)^ for large x: when this is constant it is equal to C^. We follow 
the strategy of eqs. (|5.8l5.9p . which provides more clean numerical results. 

It is expected that this procedure should work less well for b very low and a very large because in 
this case the asymptotia is reached for very large x and the numerical solution accumulates errors. 
However, for 6 — and for any a we already have the exact result derived analytically, eq. (|4.3p . 

In practice, this works well for I — 0, see for instance Fig. 1. 

For I = 1 the quantity Fi{x)'^ + Fi{x — 7r/2)^ may sometimes continue to show decreasing oscil- 
lations: here it is convenient to take into account the sub-leading term in the asymptotic expansion 
(jS.lOp . which we know to be the free wave function up to the phase shift Si. Therefore, the improved 
version of (|5.10p for Z = 1 is: 

„ / , It^ . . cos( X — ^-^ + Si) . 
Vi{^).^oo ^C-{ sin(x -- + Si)+ ^) (5.11) 

To get C^, one defines, in three steps, with Fi{x) :— (pi{x), 

m := (8V„j2)) (^i(^ + V4)-+F,(.-./4)-), 

h{x) := (tt^ - 47ra;^)(fc(a; + 7r/4) + fc(x - 7r/4), 
h{x + 7r/4) + h{x - 7r/4) 



j{x) := -i 



(7r2-16a;2)(37r2 + 16(l + x2))' 
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Figure 1: Fa{x)^ + Fo(x - 7r/2)2 for a = 1/100, = lGev,m= lO^Gev,v = 2 x 10^ 

For X large (say x > 30), j{x) quickly converges to a constant which equals C^. (For I > 1 the free 
wave function is more complicated and one should do more steps). It must be said that the results 
of the improved procedure differ little from what could be obtained simply by finding by eye the 
average of the oscillations of Fi{x)'^ + Fi{x — 7r/2)^, see Fig. 2. 
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0.0004 
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Figure 2: Fi{x)^ + Fi(x — 7r/2)^ (left) and j{x) (right) for the same parameters of Fig.l 



The case of the S / = wave has been discussed at length in the literature, and it has been found 
a resonant pattern, see in particular refs. [B] and [S]. We have verified that we get the same results. 

Here we present some numerical result for the P Z = 1 wave, which also shows a resonant pattern. 

In Fig. 3 we show the enhancement eq. (|5.7p ior I — 1, taking the values of the parameters used 
in the numerical evaluations for the S wave in ref. [5] and also reported in ref. 0, that is a = 
1/30, /i = 90Gev, as a function of m (expressed in Gev) for Vgingie particle — 10~^, 10"'*, 10~^. 

In Fig. 4 we show the enhancement eq. (|5.7p for / = 1 , taking the values of the parameters used in 
the numerical evaluations for the S wave in ref. 8J for the range of ref. [S], that is a = 1/100, fi = 
IGev, as a function of m (expressed in Gev) for Vsingie particle = 10"'^, 10""*, 10~^. 
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A The singularities in qq of the integrand at the r.h.s of 
eq.(l23]) 

Consider the iteration of eq.(2.1). At the first order we put Aq in the place of A in the r.h.s. Since 

p' is not touched in the integral equation, we can directly take its mass-shell value: 

Pg = uj{p),'p' = uj{p)n , U!{p) = \/p'^ + m? (treating the final standard-model particles as massless). 

The singularities of Aq come from the denominator of the propagator exchanged between the 
vertex x{'li)oi{p'i) s-nd the vertex x{l2)o,{p'2) (we take a massive propagator with a mass m of the 
order of -maybe equal to- m): 

{q-p'f -ni^ +ie= [q^ - {u{q,p) + uj{p)) + ie][qo + {u{q,p) ~ oj{p)) ~ ie] (A.l) 
where u{p, q) — \J u){p)'^ + Lb{qY — 2g • nuj[p) and Lj{q) = \/ q"^ + m^. 
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Therefore the integration over qq in eq. (|2.2p is (remember Pq = Lu{p), P = 0) 



I 



dqoN{qo,q,p} X — — ^ . ■ l , / / ^ , / ^^ (^'2) 

[go - {uj{q) - uj{p)) + ie\[qo + {uj{q) + uj{p)) - le] 

1 



[qo - + uj{p)) + ie][qo + {uj{q) - i^{p)) - ie] 

1 

[qo - {u{q,p) + +ie][qo + (u{q,p) - i^{p)) - ie] 



where N is some numerator, polynomial in go- We do the integration on qo in eq.(A.2) by closing 
the contour in the lower half-plane. 
The relevant poles are 

{l)qo = Lo{q)~Lo[p) (A.3) 

(2) qo = uj{q)+uj{p) 

(3) go = V^iPy +'^{qf - 2(f • nw(p) 

By computing the residue of the countour integration, it is seen that only the residue of the pole 
(1) contains in the denominator the small factor uj{q) — uj{p) — > (for m —> oo), providing the 
leading contribution in the non-relativistic limit. Therefore, in this limit it is tantamount to do the 
contour integration disregarding the singularities of Aq, that is of A at the first order iteration. 

This remains true for the next iterations since the instantaneous part of the VB propagator (that 
is the step which is making the ladder) is qo independent. 

B Derivation of the formula eg. (12. 141 ) 

Eq. p.l2p can be recast as an integral equation (remembering our choice fc^ — ie) 

r — 1 2771 CKC 
0-(f) = e*^-'^ + 77,-(f) where ^ _g2 _ p + ,g ' ^k(^ (^-1) 

Therefore 

' lim(p-^/) ^fr7^^^ ^^(p-fc) (B.2) 



(27r)-^ p-f k? — p^ — ie 



{2itY p-p k"^ — p"^ — ie 

The second term in the r.h.s. is zero because ^ '^]l_p2^l}^^ is finite (and also its integral over k 
is finite). We limit ourselves to check the first order in perturbation theory: 



= -^7 ■■ ■ = 47r^m,a / dqe^^ "- ^ (B.3) 

' ~d^~k^ + ie r J q-k^ + i({q-ky+fi'^ 

therefore 

1 J dfe~'P-%{r) _ iTT^mra 1 . . 

This is finite for generic p and k. Moreover its integral over k is also finite: by first doing the angular 
integration we get 



(27r)-^ J k'^ — p^ — ie p J -oo (^^ — — ?e)^ 
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This is seen to be finite by closing the contour in the upper plane, the singularities being a double 
pole at k = p + ie and branch cuts at A: = ±p + i^. 
Therefore there are no hidden S functions: 
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